Analytical study of four-wave mixing with large atomic coherence 



E.A. KorsunskyW, T. HalfmannW, J. P. Marangos^ 2 - 1 , and K. BergmaiW 1 ) 
W Fachbereich Physik, Universitat Kaiserslautern, D-67663 Kaiserslautern, Germany 
2 ' Physics Department, Blackett Laboratory, Imperial College, London SW7 2BZ, United Kingdom 

(Dated: February 9, 2008) 

Four-wave mixing in resonant atomic vapors based on maximum coherence induced by Stark- 
chirped rapid adiabatic passage (SCRAP) is investigated theoretically. We show the advantages of 
a coupling scheme involving maximum coherence and demonstrate how a large atomic coherence 
between a ground and an highly excited state can be prepared by SCRAP. Full analytic solutions of 
the field propagation problem taking into account pump field depletion are derived. The solutions 
are obtained with the help of an Hamiltonian approach which in the adiabatic limit permits to reduce 
the full set of Maxwell-Bloch equations to simple canonical equations of Hamiltonian mechanics for 
the field variables. It is found that the conversion efficiency reached is largely enhanced if the phase 
mismatch induced by linear refraction is compensated. A detailed analysis of the phase matching 
conditions shows, however, that the phase mismatch contribution from the Kerr effect cannot be 
compensated simultaneously with linear refraction contribution. Therefore, the conversion efficiency 
in a coupling scheme involving maximum coherence prepared by SCRAP is high, but not equal to 
unity. 



I. INTRODUCTION 

Nonlinear frequency conversion processes in atomic or 
molecular gases have attracted much attention since the 
early days of nonlinear optics. The interest is mainly mo- 
tivated by the possibility to generate coherent radiation 
in the XUV and VUV frequency range, where there are 
no transparent nonlinear crystals. However, the conver- 
sion efficiencies are usually relatively poor due to small 
nonlinear susceptibility for the generation and difficulties 
to prepare proper phase matching conditions. Approach- 
ing atomic resonances enhances the nonlinearity, but at 
the same time absorption, linear refraction and unwanted 
nonlinear phase shifts increase rapidly. 

Recently, a new technique has been put forward which 
substantially improves the nonlinear-optical properties of 
a medium. The technique, usually referred to as "non- 
linear optics with maximum coherence" , is based on the 
preparation of all atoms in the medium in the same co- 
herent superposition of two states |1) and |2) with equal 
probability amplitudes [0. From the classical point of 
view, coherently prepared atoms represent an ensemble 
of dipoles all oscillating at the same frequency with 
the same phase and with maximum amplitude. If a ra- 
diation field of frequency 0J3 is applied to the medium, it 
will beat against this strong local oscillator to produce 
the sum- or difference frequency W21 ± <^3- In this case, 
the nonlinear susceptibility of the generation process is 
large (in fact, it is resonantly enhanced) and is of the 
same order as the linear susceptibility. Therefore, com- 
plete conversion occurs within an optical length smaller 
than the coherence length. Consequently, requirements 
for the phase matching are substantially alleviated and 
the influence of density-dependent detrimental effects is 
minimized. 

In first proposals and experimental implementations 
jl], ^| , maximum coherence was established in a lambda- 
type coupling scheme with ground |1) and lower ex- 



cited 1 2) states using stimulated Raman adiabatic pas- 
sage (STIRAP) |J. However, in a lambda-type coupling 
scheme involving one-photon transitions the generated 
radiation LO21 ± ^3 cannot reach far into the vacuum- 
ultraviolet spectral region Multi-photon excitations 
may not be used in order to reach higher lying states, 
because laser-induced Stark shifts, which are intrinsic to 
multi-photon transitions, perturb the adiabatic popula- 
tion dynamics and prohibit the preparation of a maxi- 
mum coherence Q]. 

In the present paper, we investigate the use the Stark 
Chirped Rapid Adiabatic Passage (SCRAP) technique 
& B to prepare maximum coherence. In SCRAP, a 
pump laser couples a thermally populated state (most 
likely the ground state) to an excited state and a second, 
strong radiation pulse to induce a dynamic Stark shift. 
This Stark shift serves to sweep the atomic transition fre- 
quency through resonance with the pump laser frequency, 
mediating thereby an adiabatic passage of population be- 
tween two states. Provided the dynamic Stark shifts, 
induced by the second laser are larger than the shifts, 
induced by the pump field, any multi-photon transition 
may be used for the pump transition. For a two-photon 
pump transition, it is therefore possible to create coher- 
ence between a highly excited state and the ground state. 
If ultraviolet radiation is used for the pump laser, coher- 
ence between states with energies up to 10 eV may be 
efficiently created, permitting the generation of VUV ra- 
diation well below 150 nm. 

The potential of the nonlinear optics with maximum 
coherence has been demonstrated for the regime of unde- 
pleted coherence and hence for an undepleted pump field 
Q. We consider here a process of difference- frequency 
mixing involving a two-photon transition |1) — |2) res- 
onantly excited by a strong pump field with frequency 
u>\ and a dipole-allowed transition |2) — |3) excited off- 
resonance by an "idler" field with frequency u>2 (Fig. 1). 
In this case, energy for the generated field is taken only 
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from the pump field, which at the same time participates 
in the preparation of the coherence. Therefore, it will 
unavoidably be depleted if considerable conversion effi- 
ciency is expected. It is the aim of the present work to 
clarify, how the conversion proceeds when the pump field 
is depleted, what fraction of the total energy of the pump 
field may be transferred to the generated (and idler) field, 
and which parameters are needed in the specific case of 
coherence preparation by SCRAP. To this end, we solve 
the nonlinear propagation problem taking into account 
the pump field depletion. The solution of such a nonlin- 
ear problem is particularly challenging for pulses and is 
in general possible only numerically. In order to obtain 
analytical solutions we apply the so-called Hamiltonian 
approach || ||, [l(J which allows for a solution in a wide 
range of physically relevant situations. An essence of 
this formalism is to reduce a set of Maxwell propagation 
equation to canonical Hamilton equations of classical me- 
chanics, which admit several integrals of motion. Addi- 
tionally, this approach allows to analyze phase matching 
conditions taking into account intensity-dependent (Kerr 
effect) contributions, which are not present in the sim- 
ple treatment of undepleted coherence. This approach 
is especially useful under adiabatic conditions, i.e. when 
the atoms are excited by the laser pulses in such a way 
that they remain in the same instantaneous eigenstate 
of the interaction Hamiltonian during the entire process. 
This is the case in the atomic superposition prepared by 
SCRAP, as it is discussed here. 

Obviously the large atomic coherence should be main- 
tained for the duration of the conversion process. Since 
the coherent superposition includes a highly excited 
state, this requirement restricts the conversion to a 
regime with laser pulses of duration shorter than the 
natural lifetimes in the system. This is however not in 
contradiction with the adiabaticity of SCRAP. The adi- 
abatic approximation requires a slow rate of evolution as 
compared to the frequency separation of the adiabatic 
eigenstates. This results usually in a requirement for the 
product of the pulse duration and the Rabi frequency of 
the radiation field to be much larger than unity. Thus, 
for sufficiently intense fields, the process can be adiabatic 
even for short pulses. In the present paper, we consider 
the frequency conversion of short pulses and disregard 
spontaneous relaxation processes. 

The paper is organized as follows. In Sec. II we discuss 
general features and advantages of nonlinear frequency 
conversion processes with respect to atomic coherence. 
The coherence is assumed to be undepleted. In Sec. Ill 
the preparation of large atomic coherence by SCRAP 
is described. Sec. IV and the Appendix outline the 
Hamiltonian approach in nonlinear optics. Making use of 
this formalism we derive analytic solutions for frequency 
conversion processes involving maximum coherence in 
Secs.V, VI. Section V is devoted to solutions at small 
propagation distances in which the pump field is not de- 
pleted. Conclusions of Sec. II are confirmed and specific 
phase matching conditions for the SCRAP method are 



derived. In Sec. VI full analytical solutions are obtained 
taking into account depletion of the driving field. In Sec. 
VII we discuss the conditions for phase matching with 
respect to the compensation of phase mismatch, induced 
by linear refraction and the Kerr effect simultaneously. 
Finally, in Sec. VIII, we consider the spatio-temporal 
dynamics of the generated radiation pulse as well as the 
evolution of the total conversion efficiency. Conclusions 
are presented in Sec. IX. 

II. FREQUENCY CONVERSION WITH 
UNDEPLETED ATOMIC COHERENCE 



The (pulsed) e.m. field propagating in an ensem- 
ble of three- level atoms (Fig. 1) is assumed to consist 
of three components with carrier frequencies lo\ , LO2 and 

E(z,t) =^2(Sj(z,t)exp(-i(u j t-'k j r)) +c.c.) . (1) 

i 

Here |kj| = rijUjj/c with the refractive index rij at fre- 
quency uij describing refraction due to levels outside the 
three-level system in Fig. 1. The radiation pulses are 
supposed to be shorter than the relaxation times in the 
atomic system. The waves k 2 and k 3 propagate at small 
angles with respect to the vector ki (the z axis). 

1 > 




FIG. 1: Resonant four- wave mixing with maximum coherence 
between |1) and |2) adiabatically prepared by a strong drive 
uji field and a "Stark-shifting" far-off-resonance uj s field. 

In the approximation of slowly varying amplitudes and 
phases Maxwell's propagation equations read in a moving 
frame 



where £j and Vj (j = 1, 2, 3) are functions of the coordi- 
nate z and the retarded time r = t — z/c. 
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Vj are the components of the medium polarization: 



Tj exp(- 



■i(u>jt — kj-r)) + c.c. 



The pulse at frequency lu s ("ac Stark-shifting pulse") 
propagates along the z axis and is far detuned from any 
atomic state. Therefore, we assume that its intensity 
does not change along the propagation path. In reality, 
the presence of the pulse at u> s leads to generation of 
2u>\ ± lu s frequency components. However, an efficiency 
of their generation is much smaller than that for the lo 3 
field due to the off-resonant character of the interaction 
of the u> s pulse with the medium. Moreover, the presence 
of these (weak) components does not influence processes 
considered in the present paper. Therefore, we disregard 
both the change of the uj s pulse intensity and the gener- 
ation of 2cji ± oj s components. The purpose of using the 
uj s pulse is discussed later, in Sec. III.C. 

First we consider the situation of resonant nonlinear 
optics with maximum coherence for constant amplitudes 
ci and C2 of the states |1) and |2), prepared and main- 
tained by a strong pump field at w±. The back-action of 
the atoms to this field is disregarded and thus the corre- 
sponding coupling does not need to be taken into account. 
In this regime the number of photons in the preparatory 
field(s) must be much larger than the number of atoms 
in the propagation volume. Additionally, the prepared 
atomic coherence remains undepleted if the number of 
generated u> 3 photons is much smaller than the number 
of the atoms in the relevant volume. 

The components of the medium polarization can be 
expressed in terms of the atomic probability amplitudes 
c„ in levels |1), |2) and |3). After substitution into the 
one derives the field propagation 
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(3) 



(4) 



with kj (j — 2, 3) being the projections of kj on the z- 
axis. The wave vector fc 12 of the atomic coherence c*C2 
is related to the wave vector k\ of the pump field 0, 
[llj. In the case of two-photon excitation of the |1) — > 
1 2) transition, considered in the present work, we have: 
k u = 2k\. 

When deriving Eqs. (||) and ([!]), we disregarded a con- 
stant phase of the atomic transition loop, and assumed 
large detuning |A 3 | ^> Qj, with f2 2 and fi 3 being the 
Rabi frequencies for transitions |2) — |3) and |1) — |3), 
respectively: 



\ d 3 £ j\ 

2h 



(7) 



Equations (||) and (||) are linear differential equations, 
which can easily be solved. We consider the case in which 
no £ 3 field is incident on the medium, £ 3 (z = 0) = 0. 
Introducing the normalized intensity (photon flux) 



>)j 



c\£i 



I, 



and the coupling strength 



2irujjd 2 j 

he 



(8) 
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the solution of equations (|J) and (||) reads: 



m{ z ) = V20 cosh KZ\ 1 - — — 



Ak' 

2k 



1 Ak ' 

' Ak' 



■ sinh kz\ / 1 



mo 




sinh kz\ / 1 



Ak' 
~~2k 



111) 



where 7720 = V2( z — 0) is the photon flux at the entrance 
to the medium. We have introduced the conversion coef- 
ficient k: 



P12, 



(12) 



and 



where N is the density of active atoms, d 
dipole moments of transitions |2) — * |3) (|1 
is the frequency detuning indicated in Fig. 



^3 — W31 , 



«2(3) 

r 



are 
|3», 



the 
A 3 



(5) 



with uj n i denoting the transition frequencies between the 
corresponding levels. The atomic coherence between the 
states |1) and |2) is pn — | ci j , and Ak is the "residual" 
(background) phase mismatch: 



Afc = ku 



(6) 



Afc' = Afc 



N p 3 \ Cl \ 



M2 |c 2 1 



(13) 



is the total phase mismatch, including the background 
value Afc and the contributions from resonant transitions 
|1) - |3) and |2) - |3) 

The solution Eqs. (|l(j) and (^lj) shows that there 
is parametric gain (exponential growth of intensity) for 
both LO2 and ui 3 waves with the rate n if the phase mis- 
match is compensated, Afc' sa 0. Since this rate is pro- 
portional to pi 2 , it is obviously advantageous to prepare 
atoms with large coherence on the |1) — > |2) transition. 
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M3 \ci\ 



M2 \C2\ 



Phase match, 

2Afc 

~N~ ~ A 3 

can be achieved in several ways: (i) by tuning the wave 
vector fc 2 i of the atomic coherence (e.g., via the detun- 
ing A 2 of the pump field, as in Ref. pi |ll|), (ii) by 
introducing a small angle of the u>2 wave propagation di- 
rection from the z-axis, (hi) by selecting the appropriate 
detuning A3, and/or (iv) by preparation of atoms in a 
superposition with suitable amplitudes ci , c 2 . 

The resonant contributions to the phase mismatch, Eq. 
(|l3|), are usually the dominant ones over the residual 
value Afc. When phase matching is not maintained, the 
quantity 



1 - (Afc'/2/t) 



(^3 |ci| 2 ~ M2 \ci\ 



is negative, i.e. no parametric gain but periodic change 
of intensity along the propagation path takes place: 



mo 



(^3 |ciT 



N 



M2 |c 2 | 



fJ-3 \Cl\ 



Pi |C 2 | 



2A, 



(14) 



Moreover, in this regime, a substantial transfer of en- 
ergy occurs for the maximum coherence, p±2 « 1/2, case, 
whereas it is small for the regime of conventional non- 
linear optics (weak excitation, | C2 1 2 <C |ci| 2 ~ 1). The 
amount of converted energy is larger for maximum coher- 
ence than for the conventional nonlinear optics by a fac- 
tor of the order of |c 2 |~ 2 1- We stress that the assump- 
tion of undepleted coherence assumes that the number of 
generated uj 3 photons is much smaller than the number of 
photons in the preparatory ficld(s). Correspondingly, the 
total efficiency of energy conversion (from preparatory to 
generated helds) is very small in this regime. 

The regime of undepleted atomic coherence corre- 
sponds to the classical picture of frequency mixing in 
which the atoms play the role of a local oscillator (fre- 
quency UJ21) with the "probe" L02 field beating against 
it to produce the difference (or sum-) frequency uj^ — 
UJ21 i 0J2- Such a process is obviously more efficient for 
a strong local oscillator, i.e., for large pi2- Thus, the 
preparation of a large atomic coherence is favorable for 
frequency conversion in atomic gases. 



III. PREPARATION OF MAXIMUM 
COHERENCE BY SCRAP 

A. Atomic parameters 

First we discuss the specific parameters to be consid- 
ered in the coupling scheme discussed here (Fig . 



The Rabi frequencies of single-photon transitions Qj 
(j = 2, 3) are related to the photon flux ly, Eq. (||), via 
the coefficients pj, Eq. (^), as 

{ -'-, = \ Wh- 

The phase of the Rabi frequency f2i for a two-photon 
transition |1) — > |2) is equal to 2tp\, and the module is 
proportional to the intensity: 



1 2 



(15) 



where pi is the transition coupling constant, and 
tt„„i (luj) is the matrix element of an atomic polarizabil- 
ity tensor: 



TlCtr, 



Pi = 



c 



/// 



ai2 (wi) , 

(n\ d \m) (m\ d \n') 



(n\ d \m) (m\ d \n'} 



(16) 
(17) 



(w m — LO n ' 1 ) + LOj 



with huj n i being the energies of the resonant states \n) 
(n = 2,3, with un = 0), hio m the energies of the (vir- 
tual) states I to), and (n\ d \ m) the dipole moment matrix 
elements for transitions \n) — > |m) (n — 1,2,3). 

The frequency detunings A„ (n — 2,3) include the 
"static" detuning <5„o, a c Stark shifts n jVj induced by 
the ujj (j — 1,2,3) fields, and the shifts S n — flnsVs in- 
duced by an intense far-off-resonant " SCRAP laser pulse" 
with frequency lo s and photon flux r] s : 



J = l,2,3 



/3 njVj , (n = 2,3), (18) 

(19) 



$30 
$20 



2ft LUj 



(a nn (ujj) - an (wj)) 



(20) 



w 3 — w 31i 
2wi - UJ 2 1- 



It is important to note some essential relationships be- 
tween the atomic parameters used in the present work. 

For atomic media the off-resonant (background) con- 
tributions to the refractive index rij are expected to be 
of the order of [HI: 



1 + 2nNan (ujj) . 



(21) 



The residual phase mismatch Afc = 2ki — k 2 — k^ is there- 
fore of the order of 



Afc : 



2nN 



{2u-ia.ii (ojx) - lu 2 uii (w 2 ) - W3CH1 (^3)) ■ 

(22) 
that the 



It follows then from Eqs. (|16|), (|17]), (gOJ), (£j 
quantities: 

[ii ~ [3 nj ~ Ak/N 
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are all of the same order of the ma gnit ude. 
Further, we have from Eqs. (g), @, @: 
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— \(u> m - UJnl) ±UJj\ > S 30 , 520, \0,j\ (23) 

Ml 

The last inequality is implied by the resonant three-level 
model of the atom. The validity of this inequality justifies 
the use of the rotating wave approximation. 

Finally, we present values of the above constants for 
a real atomic scheme, which can be used to drive the 
generation of short wavelength radiation. We consider 
a coupling scheme in Kr with the states: |1) = 4p 6 1 S* 
(ground state), |2) = Ap 5 hp [0, 1/2] (94 093.7 cm,- 1 ) and 
1 3) = 4p 5 5s [1,1/2] (80 917.6 cm" 1 ). The two-photon 
transition between the ground and the excited state in 
Kr is known as an efficient transition for conventional 
four- wave mixing schemes. The scheme discussed here 
has e.g. been used in experiments on VUV generation 
assisted by electromagnetically induced transparency |ll| 

H . The pump field at 212.55 nm excites the two-photon 
transition |1) — |2), the idler field at 759 nm is tuned near 
the single-photon resonance of the transition |2) — |3), 
and the field generated on the |3) — |1) transition has a 
wavelength of 123.6 nm. The coupling strength of the 
single-photon transitions used in this scheme is : [12 = 
3.507 x 10 -2 cm 2 x s _1 and ^ 3 = 0.441 x 10~ 2 cm 2 x s -1 . 
The coupling strength /ii of the two-photon transition 
and the ac Stark coefficients (5 2 j can be estimated as \i\ ~ 

I. 8 x lO" 16 cm 2 , fci ~ 3.7 x lO" 17 cm 2 , /3 22 « 2.2 x 

,2 



-17 



cm 2 , ^23 ~ 6.4 x 10 17 cm 2 . The residual phase 
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mismatch for this scheme have been measured 
value is: Ak/N = 4.8 x 10" 



" 17 cm 2 . 



B. Interaction Hamiltonian, eigenvalue equation 

In rotating-wave approximation, the light-atom inter- 
action Hamiltonian is given by: 

H = -fc[A 2 |2)(2|+A 3 |3)(3|] (24) 

-mi |i) (2| + nn 2 e i<p \2) (3| + nn 3 \i) (3| + h.c, 



where the multiphoton resonance condition 



UJ3 = 2oj\ 



L02 



(25) 



has been used. 

In the present paper, we consider adiabatic light-atom 
interaction processes, i.e. the atomic system can be as- 
sumed to follow the evolution of the instantaneous eigen- 
states. If, for example, the atomic system is at some 
initial time to in the nondegenerate eigenstate IV'o(^o)) of 
the interaction Hamiltonian, i.e. 



(26) 



(which is usually the ground state of the atoms), it will 
remain in this state jV'o) at ah times. Eq. (p6[) yields the 



characteristic equation for the eigenvalues: 

Ac (A 2 + A ) (A 3 + A ) - (n 2 + n 2 + n 2 ) a 

-ttf A 3 - O3A2 = -20^2^3 cos <p, (27) 
where the relative phase ip of the elm. waves is: 

ip = 2<pi — ip2 — <f3 — Akz, (28) 
which includes the residual phase mismatch Afc. 



C. Preparation of maximum coherence by SCRAP 
procedure 

In what follows, we will concentrate on a regime of 
nonlinear optics with large coherence between the states 
|1) and 1 2). To prepare such a coherence efficiently, we 
suggest to use the Stark Chirped Rapid Adiabatic Pas- 
sage (SCRAP) method [|[ g{. In the particular example 
of three-level system in Fig. 1, this procedure can be 
realized when the idler L02 field is far detuned from the 
resonance with transition |2) — |3) (i.e., the static detun- 
ing (530 is much larger than other parameters including all 
the Rabi frequencies and detuning A2). In this case, the 
adiabatic (dressed) state that asymptotically connects to 
|1) for t — > —00 is given by: 



|1> 



An 



the with corresponding energy eigenvalue 



A w -^A 2 



1 



4f! 2 . 



|2) , (29) 



(30) 



As it is obvious from these formulae, population can 
be prepared from the bare state |1) in the dressed state 
IV'o)) if at the beginning of the interaction A2 — * +00. 
The state |^o) will project completely onto the target 
bare state |2) at the end of the interaction, if A2 — > —00. 
Thus all the population can be transferred from the 
ground to the excited state via the dressed state. For 
pulses with duration in the nanosecond range, this trans- 
fer process can be implemented experimentally by sweep- 
ing the atomic transition frequency with an additional 
laser pulse of frequency tu s (see Fig. 1) inducing dynamic 
Stark shifts (SCRAP). The pump and Stark shifting laser 
pulses have to be delayed. Otherwise the effect of rapid 
adiabatic passage will occur twice, once in the rising, sec- 
ond in the falling wing of the Stark shifting laser. 

Fig. H shows the population and coherence dynamics 
induced in an atomic system driven by laser pulses in 
SCRAP configuration. As described above, population 
is transferred completely from the ground to the excited 
state, as the atomic transition frequency is swept in the 
falling edge of the Stark shifting laser pulse through res- 
onance with the pump laser (middle frame). The coher- 
ence P12 induced in the system during the interaction 
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While the coherence, induced by SCRAP, is not per- 
manent, a slight modification of the process permits the 
preparation of a long lasting maximum coherence (so 
called half-SCRAP) g. In this configuration the pump 
laser is tuned to resonance, i.e. the static detuning 
^20 = 0- Fig- |H shows the population dynamics and 
the coherence induced in this case. At earlier times we 
have A2 3> r2i so that the state |-0o) coincides with the 
ground state |1). When fii increases, the adiabatic state 
\ipo) evolves in a superposition |1) and |2). At the end 
of the interaction the situation Sli 3> IA2I is reached, 
and the adiabatic state corresponds to "maximum co- 
herence": \ip Q ) = (|1>-|2»/V2. This regime requires 
to fix the static detuning S20 to zero with sufficient ac- 
curacy, but it establishes the enduring coherence needed 
for phase matching. 



FIG. 2: SCRAP procedure producing complete population 
transfer and a pulse of large atomic coherence. Time evolution 
of the detuning A2 and the two-photon Rabi frequency fii 
(top frame), the populations of the |1) and |2) states (middle 
frame), and the coherence pi2 (bottom frame). 



reaches a maximum of 1/2 when the population is dis- 
tributed equally between the bare states |1) and \2). The 
system is prepared in maximum coherence. This hap- 
pens, however, only at one instant of time, which results 
in a pulse of large atomic coherence p±2 (t)- We note 
that the transient large coherence p\ 2 occurs also when 
the pump and Stark shifting laser pulses coincide. As we 
show later, however, such regime leads to quite low overal 
conversion efficiency due to phase matching reasons. 



IV. HAMILTONIAN APPROACH FORMALISM 

In this work, we use an approach which does not 
require explicit expressions for the atomic amplitudes 
H H |fO|].The main advantage of this approach is the re- 
duction of Maxwell propagation equations (§) to the form 
of canonical Hamilton equations of classical mechanics 
involving action and angle variables J and ip (see Ap- 
pendix). Here ip is the relative phase of e.m. waves, Eq. 
(p8|), and the variable J(z) characterizes the amount of 
energy exchange between the waves and has the initial 
condition J(z — 0) = 0: 



Vi( z ) = Via - 2J( Z ), 
m{z) = V20 + J{z), 

m{ z ) = V3Q + J(z). 



(31) 
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After some algebra (see Appendix), the Hamilton 
equations can be further reduced to yield an implicit so- 
lution for J{z): 



S(J') 



dJ' 



(32) 



where both functions R (J) and S (J) are polynomials in 
J. 

For the generation of the 013 mode from vacuum: 7730 = 
0, the functions R (J) and S (J) take a form: 

R = ^lp 2 p, 3 J(r, w -J) 2 { m + J) (33) 

-(A 1 + A 2 J + A 3 J 2 ) 2 J 2 , 
S = a + a\J + a 2 J 2 . (34) 



FIG. 3: "Half-SCRAP" procedure leading to a permanent co- 
herence after the interaction. Time evolution of the detuning 
A2 and the two-photon Rabi frequency fii (top frame), the 
populations of the |1) and |2) states (middle frame), and the 
coherence P12 (bottom frame). 



As shown in the Appendix, cf. Eq. ( A21 ), the coefficients 
A m and a m describe the linear and nonlinear refraction 
coefficients of the medium. 

Equation (|3^) matches a one-dimensional finite mo- 
tion of a pendulum in an external potential. The allowed 
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range of J, corresponding to the region of classically al- 
lowed motion of the pendulum, lies between zero and the 
smallest positive root J\ of the polynomial equation: 



R(J) = 0. 



(35) 



The second term in expression ( p3| ) for R(J) is never 
positive, so the smallest positive root of the polynomial 
is bounded by r/io- This reflects the fact that the conver- 
sion process stops when the energy of the pump field is 
entirely depleted. In order to reach this limit and thus to 
attain maximum conversion efficiency, the second term 
in ([53j) should be small, which corresponds to negligible 
phase mismatch : 



Ai +A 2 J + A 3 J 2 



0. 



In order to see which values are required to approxi- 
mately satisfy this condition, we have to analyze the co- 
efficients A m . 

The coherence preparation process requires large static 
detuning S 3 q and small ac Stark shift induced by the 
idler 0J2 wave: LI2V20/830 ^lVw (see discussion in Sec. 
III.C). Taking into account the relation /X2/M1 ^ ^30; Eq. 
( p3| ) , the latter requirement restricts the intensity of the 
W2 wave: 7720 "C rj 10 . Since in the down-conversion pro- 
cess considered here the energy is taken only from the wi 
wave, this condition does not impose a real limitation. 

Under these conditions, the non-vanishing coefficients 
A m and a m are given by 



A, 
A 2 

Oi 
where 



-qS 30 (2A + 6 2 + faiVw) 
-/i 2 A - ^3 (A + 6 2 + Auf7io) , 
q 2 S 30 + H2Q + M3 (q ~ /?22 - #23 
6 30 (2A + 62 +/?2i?7io) , 

- (M2 + M3) , 



q = 2Ak/N 



(36) 

2/?2i) , (37) 
(38) 
(39) 



(40) 



The eigenvalue A is a con stant of motion which can thus 
be found from Eq. (A20): 



A = A = -- (62 + /3 2 i?7io) 



+ y(S2+(32iViof +^Ww (41) 



SOLUTIONS FOR UNDEPLETED PUMP 
FIELD 



V. 



We first consider the solution of Eq. (|32j) for small 
density-length products Nz in the case of low conversion 
efficiencies, i.e. the intensity of the generated field is 
assumed to be much smaller than the intensity of the 
pump u>i wave: 



773 = J < ?7io- 



We can then neglect the term A2J in expression ( |32 
and the term a\J in expression (|3~i|). In this case the 
solution of Eq. (B3) has a simple form: 



»/20 



' Afc' 
^ 2k 



sinh 




Afc' 
"2k" 



Afc' 



= Afc 



(530 2A + ^2 + /?2i?7io ' 
2 a 

N Li 3 (A + (5 2 + IhiVw) + M2A 



(5 30 (2A + (5 2 + /3 2 i?7io) 



(42) 
(43) 

(44) 



which is similar to Eq. (|llj) obtained under the assump- 
tion of constant probability amplitudes of the states |1) 
and 1 2) (cf. Sect. II). 

The optimum conversion (parametric gain with large 
rate k) occurs when the phase mismatch Afc' is compen- 
sated. Taking into account Eq. (^) for A, the condition 
for phase matching reads: 



N Ak = 



1*3 + M2 

2(530 

M3 - M2 



((5 2 + foiVw) 



2(5 



3(1 



\J (82 + An^io) 



(45) 



It is very important to recognize that the parameters (5 2 
and 7710 are time-dependent (pulsed) in the SCRAP pro- 
cess. Therefore, the r.h.s. of Eq. pq ) is time-dependent. 
Hence, it is impossible to phase-match the generated and 
the pump waves for the duration of all stages of the light- 
atom interaction process. This fact has a detrimental 
effect for the full SCRAP procedure where the pulse of 
large coherence is produced (see discussion in Sec. III.C 
and Fig. |^). However, in the half-SCRAP case with the 
permanent large coherence (Fig. ^), there are two rela- 
tively long time intervals, in which the phase matching 
condition does not depend on time. 

At the early stage of the SCRAP process, when 62 3> 
Mi^iOj the condition is: 



N 



-Afc : 



M3 

ho ' 



(46) 



and the nonlinear conversion coefficient takes the form: 

Mi ??io 



Kl « k - 



(47) 



with the "maximum coherence" conversion coefficient kq: 



N V^3 



''30 



(48) 



For later times, when the Rabi frequency of the pump 
field exceeds the detuning: jXirjiQ 3> 82 and the adiabatic 
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state corresponds to the maximum coherence superposi- 
tion, the phase matching condition becomes: 



— Ak - ^ 3 +M2 _ M3 - fin 

N ~ ~ 2S 3Q 2<5 30 y/p^ + 4tf ' 

and the conversion coefficient is 



K 2 ~ KO 



(49) 



(50) 



Phase matching according to Eqs. (1461), fl49|) can be 
performed by controlling the background mismatch Ak 
through a suitable choice of the small angle of the w 2 wave 
propagation direction from the z-axis. For radiation in 
the visible spectral range, detuning ^30 ^100 GHz and 
atom densities N ~ 10 13 4- 10 14 cm~ 3 , the value of Ak 
necessary to com pen sate the resonance refraction con- 
tributions, Eqs. (|46|), (pL9|) , corresponds to an angle of 
0.1 — 1 rar ad- 
it is obvious from Eqs. j47|), ( p0| ) that k 2 » K\. There- 
fore, it is advantageous to drive the frequency conversion 
process when a large coherence is established, thus to 
apply the pulse at u>2 when fiirjio 3> 82 and to choose 
parameters satisfying the condition Eq. (|49|). 

When phase matching is not maintained, the quantity 
1— (Ak' /2k) 2 is always negative (Ak is much smaller than 
the resonant contributions to Ak'), and there is no expo- 
nential growth but sinusoidal oscillations of the generated 
intensity with respect to Nz. However, for /117710 ^> 82 
(maximum coherence) the quantity 1 — (Ak' /2k) 2 in 

the denominator of Eq. (^2|) is of the order of unity, 
whereas for S 2 A*i??io (atoms are in the ground state) 
we have 1 - (Ak' /2k) 2 ~ (Ak' /2k) 2 > 1, and corre- 
spondingly, conversion is tiny. 

These considerations, similar to those of Sect. II treat- 
ing the nonlinear conversion with fixed probability ampli- 
tudes, demonstrate once again that large atomic coher- 
ence is preferable whatever the method of the preparation 
might be. 

Thus, in the case of good phase matching, we observe 
exponential growth of the generated intensity. The ques- 
tion arises up to which values the generated intensity 
will grow and what the limiting factors are? Since the 
pump wave will be depleted, one may also expect that 
the preparation of large atomic coherence will not be ef- 
ficient anymore. By then it is not clear how this will 
influence the frequency conversion process. In order to 
answer these questions we need to solve the complete 
propagation problem taking into account the depletion 
of the pump field. 



VI. 



GENERAL SOLUTIONS FOR RESONANT 
FOUR- WAVE MIXING 



The solution of the propagation equation (32) is deter- 
mined by the roots of cubic equation (Bq), in particular, 



by their signs and relation between their modules. Un- 
der the condition ^ 2 o <^ ^10 j the roots of ([35]) can be well 
approximated by: 

1-h 



X2 = 



X3 



I + 62' 
I-62' 

ma 1 



mo i-bj 



(51) 
(52) 
(53) 



where 



The quantities 61, 62: 
bi 



2 Ml??10v / M2M3 ' 

A 2 



2^1^/7^27*3 

determine the phase mismatch induced by linear refrac- 
tion and Kerr effect, respectively. 

Since 1720 / ?7io 1 we have in most relevant cases: 
\xa\ < \xi\ , \x 2 \. 

Evaluation of the integral in Eq. (pi) gives the follow- 
ing general dependence for x(z) = J(z) /rjio in implicit 
form: 

±k'z + xo = F[y (x),p) 



a 



r{F[ 1 (x),p}-dU[ 1 (x),n,p}}, (54) 



where xo is an integration constant, and F(j,p) and 
11(7, n,p) are the elliptic integrals of the first and third 



kind, respectively ]lf 
coefficient defined as 



is the nonlinear conversion 



a V V \xi\J s 



(55) 



The parameters of the elliptic integrals 7 (x) ,n,pas well 
as the factors r, s,d depend on the signs of expressions 
(l-&i) and (l-6§). 

Due to the condition 7720/^10 "C 1, the expression Q) 
can be inverted to give the explicit solutions presented 
below. 



A. Compensation of both linear refraction and 
Kerr effect : b\ < 1 and b\ < 1 

For b\ < 1 and b\ < 1 the relevant parameters are: 



7 (x) — arcsin - 



l x(xi + \x 3 \) 
xi (x+ \x 3 \)' 



= / *i(*a + l*a|) _gi (56) 
y x 2 (xi + \x 3 \) X\ + \x 3 \ 

d = 1, s = 1, r = \x 3 \ . 
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In this case, the solution is as follows: 
Xi \x 3 \ sn 2 [k'z;p] 



x(z) 



xicn^ [k z;p\ 



(57) 



where sn[/tz;p] and cn[rez;p] are the Jacobi elliptic sine 
and cosine functions, respectively pqj . 

In this case, the parameter p is close to unity, so that 
sn[Kz;p] — > tanh (kz) and cn[/tz;p] — s-sech(Kz). Thus, 
for small density-length products Nz such that k'z <C 
In [x\j | ac3 j ) the solution ( |57| ) is reduced to 

x(z) = x 3 sinh 2 (k'z) , 

which coincides exactly with the solution obtained under 
the condition of undepleted pump field, see Eq. (E3). 

For larger Nz, Eq. ( |57| ) has to be used. The form of 
this solution is shown by the solid line in Fig. ||. The 
maximum value of x(z) attainable in this regime is given 
by x\ corresponding to the intensity of the generated uj 3 
wave given by : 



Jrnax — ^?10 



1 + 62' 



which is of the order of t\\q. Thus, almost complete 
conversion can be achieved in this regime. This maxi- 
mum value is reached at the distance z — («/) 1 K (p) 
with K (p) being a complete elliptic integral of the 
first kind, which can be approximated for p « 1 — 
\x 3 \ (x 2 - xi) j (x\x 2 ) — ► 1 as 



A-(p)w(l/2)ln 



\§X\Xl 



(x 2 - xi) \x 3 \ 




B. Compensation of linear refraction b\ < 1, but 
large Kerr-induced refraction 62 > 1 

In this case, the parameters are as follows: 



7 0*0 = 



arcsm - 



> \X2\ (xi - x) 

x\ (x + \x 2 \) ' 



xi (\x 2 \ - \x 3 \) 

~ \l \x 2 \(xi + \x 3 \y 

= 1 — n, s = 1 



XX 

\x 2 \ 



(58) 



aiTyio 



aov — n 



^2 



The solution reads 



Xi IX3I sn 2 [k'z;p] 



^ 3 |+a;i^q^cn 2 [ K 'z;p] 



(59) 



The form of the solution is similar to the previous case, 
except for the prefactor at the cn 2 [k/z;j>] function in the 
denominator. The spatial evolution of the generated in- 
tensity in this case is plotted as a dotted line in Fig. ||. 

We observe a parametric gain at the initial stage of the 
propagation, at k'z < In ( jf^^f; 



-|Z2 



x(z) = x 3 



Xl + \x 2 \ 
2\x 2 \ 



sinh 2 (k'z) . 



(60) 



The maximum of x{z) from Eq. ( 59 ) is again given by x\ . 
However, for large "Kerr coefficient" b 2 ^ 1, the value of 
xi corresponds to an intensity of the generated ui 3 wave 
that is much smaller than 7710: 



I ) m ° 

\ 13/max l 
^2 



(61) 



It is also important to note that the conversion coefficient 
k' here is smaller than in the "compensated case" by 
a factor of s ~ (^2,3/^1^30 V~ n ) ^> 1- Therefore, the 
conversion proceeds much slower, see Fig. 0. 



C. No compensation of linear refraction : b\ > 1 

For b\ > 1 and b\ < 1 the elliptic integral parameters 
are: 



FIG. 4: Spatial evolution of J (z) at given retarded time r 
in the case of compensated linear refraction. Parameters: 
ma / tim) = 0.01, bi = 0.1. The solid line corresponds to Eq. 
©: b 2 = 0.5. The dotted line corresponds to Eq. (!§): 
62 = 8, s = 5. 



7(2) 

P 
d 



arcsm - 



l x (x 3 + \xi\) 
x 3 (x+ \x-lI) 1 



lx 3 (X 2 + \X!\) 



x 3 



y x 2 (x 3 + \xi\)' x 3 

1, s = I, r=\x 1 \. 



(62) 
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For bf>l and b\ < 1 



7(x) 



x 3 (x + |a; 2 |) ' 



P = 
d = 



1x3 I F2 



a* 



V N2I (£3 + l^i I) ' 
1 — n, s = 1, r 



1*2 I 



(63) 



In both cases p< 1 and n< 1. This permits a reduc- 
tion of the solution to the form: 



x (z) = \x^ \ sin 2 (k'z) . 



(64) 



The maximum of x (z) is \xs\, i.e. J m ax = V20/ 1 1 — | - 
Therefore, this solution demonstrates the crucial influ- 
ence of the phase mismatch induced by linear refraction. 
If this contribution to the phase mismatch is not compen- 
sated, the maximum intensity of the generated LO3 wave 
is always limited by the input intensity ^20 of the idler 
wave. 



VII. COMPENSATION OF PHASE MISMATCH 

As we have shown in the previous section, it is crucial 
to compensate the mismatch induced by linear refraction 
b\ < 1 in order to get large conversion. Only then expo- 
nential gain occurs at the initial stage of the process and 
the maximum generated intensity will be much larger 
than the input intensity 7720 of the idler 0J2 wave. At 
the same time, it is desirable to make the Kerr-induced 
mismatch as small as possible. 



A. Compensation of the phase mismatch induced 
by linear refraction 



In general, the condition b\ < 1 yields 



2/o 



< y < yo 



s/l + dl V 1 + dl 

where the "phase-matching-tuning parameter" y is: 

q&30 



y 

The value of yo: 

1 — to di 



VM2M3 



(65) 



(66) 



1 + m 



yo = 



with notations 

M2 

to = — . 

, Ihl , &2 



2/j.x 2/xi7jio' 



(y < 0) , (67) 

(68) 
(69) 



corresponds to the condition given by Eq. (|4^) where 
b\ = (complete compensation of the linear refraction) . 
At the limits 



j/i,2 = yo t 1/v 1 



d\ 



(70) 



of the desirable range of y, we have b\ — > ±1. From Eqs. 
(pl|)-(|53|) we see that the root x\ determining the max- 
imum conversion efficiency becomes very small at these 
values of y. Therefore, it is not favorable to set the work- 
ing point close to 7/1,2- 
_We stress an important consequence of the inequality 
For any 62 and /ti^io, the quantity y < as well 
as its absolute value is of the order of one: \y \ ~ 1 in the 
range given by Eq. (pq). Therefore, q ~ ^2,3/^30 and 



\l>2 



M2,3 



» 1. 



:!(] 



The Kerr-induced phase mismatch is large in the range of 
parameters (q and 630) while linear refraction is compen- 
sated in this regime. Thus it seems impossible to simul- 
taneously compensate both contributions of the phase 
mismatch. 



B. Compensation of the phase mismatch induced 
by the Kerr effect 



Condition b\<\ yields 



1 + to 



y\y + 



< 



2/ll<5; 



311 



(71) 



Due to 2//i$3o/-y/M2M3 ^ 1; the above condition can be 
fulfilled by 



\y < 



2/11^30 



M2 + M3 

i.e., by y f=a 7/3 = (or equivalently, by q « 0), and by 



(72) 



y- 



1 + m 



< 



2/ik5. 



3(1 



VM2M3 



(73) 



that is by y ~ 7/4 = — (1 + to) / \fm [or equivalently, 
by 9^30 ~ — (a*2 + ^3)1 see Eq. (E9I)] in a very small 
range ±2/ii<53o/ {fi2 + ^3)- The values for 7/3,4 are fixed 
by atomic parameters and cannot be tuned. 
Further, it is easy to show that 



J/4 < 2/1 < %J2 < 0. 



(74) 



We see that the only possibility to satisfy both b\ < 1 
and b\ < 1 is to make the value y\ close to y m and to 
tune y to the vicinity of 7/1 (and y m ). However, as we have 
discussed in the previous subsection, it is not favorable 
to work with y close to y\ since the root X\, which de- 
termines the maximum conversion, becomes very small. 
Moreover, the condition 7/1 = y m reduces to 



(1 — to) 



m 



(75) 
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what can be realized only at one instant of time, because 
d,2 is a time-dependent function, see Eq. (|69|). Then, 
only a very narrow part of the generated pulse may, in 
principle, be phase-matched to the pump waves, and any 
small deviation from the above condition given by Eq. 
((75|) will destroy the phase matching. 



VIII. SPATIO-TEMPORAL EVOLUTION. 
TOTAL CONVERSION EFFICIENCY 

As we have seen before, there are different phase 
matching regimes at different stages of the SCRAP pro- 
cess. If we choose to compensate the linear refraction 
by satisfying the relationship (|49|), then at the beginning 
when 82 3> Mi?7io the linear refraction is large: \b\\ > 1. 
When the pump intensity reaches the values /ii?7io 3> 82 
the mismatch is completely compensated. Therefore, it 
is favorable to start the conversion process, i.e., to apply 
L02 pulse at time instants when /ii?7io S> 82. However, 
the delay of u>2 pulse should not be too large since the 
conversion process is also determined by the overlap be- 
tween the pump u)\ and the idler u>2 pulses. In order to 
illustrate these processes, we show graphical representa- 
tion of our analytical results in Figs. ||, ||, 0. Figs. || and 
^ demonstrate the evolution of the generated intensity in 
the case of the permanent coherence preparation by the 
half-SCRAP, and Fig. - for the case of pulsed large 
atomic coherence induced during the population transfer 
in the full SCRAP process. 

In Fig. |^, with the idler 102 pulse arriving before 
the pump lj\ pulse, sinusoidal oscillations occure along 
the propagation path for the early part of the generated 
pulse. This corresponds to Eq. ( |64| ) for the case of un- 
compensated phase mismatch. As the pump intensity 
?7io increases and the interaction parameters get closer 
to the phase matching condition Eq. ( |49|) and \b\\ be- 
comes sufficiently small: \b\\ < 1, Eq. (p9[) is applied. 
During this time interval intensities J much larger than 
1J20 do occur. We recall that 7720 is the maximum inten- 
sity that can be obtained without elimination of linear 
refraction. However, since the Kerr-induced mismatch is 
large, the maximum generated intensity 773 = J is still 
much smaller than 7710 (it is given by 7710/^2, Eq. (|3l|)), 
and the rate of intensity growth is quite small. 

Fig. ^ shows the variation of the intensity J with r 
and z, when the phase matching takes place over the en- 
tire duration of the generated pulse. We see from Figs. 

[| that the maximum of the generated pulse always co- 
incides with the maximum of the pump lo\ pulse and not 
with that of the idler L02 pulse. This fact directly follows 
from the physics of the down-conversion process in which 
the ll>2 field serves simply as a seed wave, while the en- 
ergy is taken only from the pump u)\ field. However, the 
temporal overlap between the pump and the idler pulses 
also influences the conversion process. Initially, there is 
exponential growth with x(z) ~ X3 sinh 2 (k'z) ~ 7720/7710, 
Eq. (Bet). Thus the generated intensity is determined 



2oo ^ms 



_ *■ 



t/T, 



FIG. 5: Spatio-temporal evolution of generated intensity 
J (z,t) (normalized to the maximum of 7710) for the half- 
SCRAP preparation and linear refraction compensation at 
maximum coherence according to Eq. (fli^). The retarded 
time r is in units of the duration Ti of the pump u)i pulse, 
and the propagation distance z is in units of conversion length 
/ty" 1 for the ideal maximum coherence case, Eq. (jig]). The 
parameters are ^2/^3 = 7.95, /?2i/2/iti = 0.1 (Kr atoms), 
[12/(2^1830) = 20. The temporal profile of the Stark, pump 
and idler pulse is Gaussian with the following parameters: 
maximum of detuning and pump Rabi frequency: 8™ /Q^o = 
2, static detuning £20 = 0, rf^/rf^ — 0.005, center of the Stark 
pulse t a /T\ — —1.5, duration of the Stark pulse T 3 /Ti — 1, 
duration of the 0J2 pulse T2/T1 — 0.5, delay of the u>2 pulse 
t 2 /Ti = -1. 



mainly by the idler field: 773 ~ 7720 sinh 2 (k'z). Later, the 
maximum of the generated pulse is shifted towards the 
maximum of the uj\ pulse. This dynamics leads, in gen- 
eral, to a temporal modulation of the generated pulse. 
The best conditions are obtained when the maxima of 
the pump u>i and the idler 0J2 pulses coincide (Fig. |^). 
In this case, conversion proceeds more or less homoge- 
neously and the temporal shape of the generated pulse is 
smooth at all propagation distances. 

When the large atomic coherence persists only in tran- 
sient during the full SCRAP process, efficient generation 
occurs during a short time interval around the maximum 
of the pump pulse (see Fig. R). This interval is deter- 
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FIG. 6: Evolution of J (z, t) for the half-SCRAP preparation 
and compensation of linear refraction at maximum coherence 
by Eq. (^). Parameters are the same as in Fig. |E| except the 
delay of the U2 pulse tilT\ = 0. For better visibility, the time 
and length axes have been reversed as compared to Fig. ^. 

mined by the range where the phase matching, | £»x | < 1 , 
occurs. 

As we see, the temporal shape of the generated pulse 
is in general quite complicated. That is the output pulse 
is not transform-limited. An important quantity to char- 
acterize the conversion process is the total energy con- 
version efficiency, defined as 

/ dtuj 3 rj 3 (z,t) 
W{z) = ^- . ( 7 6) 

/ diwi?7io(2:,i) 

The evolution of W{z) is shown in Fig. || for three dif- 
ferent delays of the idler pulse in the case of permanent 
coherence. The largest conversion efficiency is obtained 
for coinciding lo\ and LU2 pulses. Moreover, the maxi- 
mum of W(z) occurs at propagation distances smaller 
than that for the case of delayed pulses. However, the 
total conversion efficiency is not substantially different 
for different delays because the maximum of W(z) is de- 
termined mainly by rf{^/b2 , i.e., by the parameters of 
the pump pulse. The thin solid line in Fig. displays 



FIG. 7: Evolution of J (z, r) for the temporally large coher- 
ence induced during the full SCRAP process. Linear refrac- 
tion is compensated for maximum coherence by Eq. (p9[). 
Parameters are: 5? = 10, 5 2 o/ n iQ = ~ 5 > t a /Ti = -1.7, 
T 3 /T\ = 2, tilT\ = 0. Other parameters are the same as in 
Fig. H. For better visibility, the time and length axes have 
been reversed as compared to Fig. 0. 



the conversion efficiency in the case of pulsed large coher- 
ence. As expected, the efficiency is much smaller than in 
the "permanent P12" case because of two reasons. First, 
the nonlinear conversion coefficient k, Eq. ([T2]), is large 
in the transient regime during a short time slot. Second, 
more important, phase matching, as discussed above, can 
be achieved only in an even shorter time interval. 



IX. CONCLUSIONS 

We have discussed the analytic solutions of a four- 
wave mixing process involving preparation of an atomic 
system driven to maximum coherence by the technique 
of Stark chirped rapid adiabatic passage (SCRAP). The 
maximum coherence permits conversion efficiencies, ex- 
ceeding the case of conventional nonlinear optics by a 
large factor. However, the conversion efficiency does not 
reach unity, because of phase mismatch due to the lin- 
ear and the intensity-dependent index of refraction. It is 
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APPENDIX A: HAMILTONIAN APPROACH 

Here we present an outline of the Hamiltonian ap- 
proach in nonlinear optics & |l This approach is 
based on the representation of the medium polarization 
P as a partial derivative of the time-averaged free energy 
density of a dielectric with respect to the electric field 
strength E p): 



P = -( N 



r dH 
dE 



(Al) 



where (...) denotes quantum-mechanical averaging, and 
H is the single-atom interaction Hamiltonian. With the 
field given by Eq. ([!]), we write: 



FIG. 8: Spatial dependence of the total conversion efficiency 
W, Eq. (f7q) , for different delays of the L02 pulse for the half- 
SCRAP coherence preparation: solid line - tz/Ti = 0, dotted 
line - tz/Ti = +1, dashed line - £2/21 = — 1. The values for 
all other parameters are the same as in Fig. |^. The lower 
thin solid line corresponds to the transient large coherence 
induced during the full SCRAP process, Fig. ^. 



P = ~ ( N Y1 Q£* exp(-i(Ujt - kjz)) + c.c. 
\ 3 i 

thus the propagation equation (0) becomes: 



(A2) 



dz 



-i2n^-N < 



dH 
d£* 



(A3) 



practically impossible to compensate both parts simulta- 
neously. Thus, the conversion efficiency gets maximum 
when the linear part of the phase mismatch is reduced 
to zero. This can be done by the controlling the resid- 
ual phase mismatch Afc through the buffer gas or non- 
collinear pulse propagation. Unfortunately, the Kerr- 
induced mismatch is large and limits the rate and the 
maximum achievable efficiency of the conversion. Still 
the maximum atomic coherence, prepared by SCRAP, 
permits efficient generation of strong short-wavelength 
radiation. In principle, the generated radiation is broadly 
tunable - the detuning 630 may be changed provided it 
is still larger than Rabi frequencies and ac Stark shifts. 
However, the phase matching condition Eq. ( ff9[) has to 
be satisfied. Therefore, tuning of £30 requires modified 
compensation of the residual phase mismatch Afc. 
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When the atomic system adiabatically follows the instan- 
taneous eigenstate \ipo), we find 



dH 
~d£* 



(W m) = Rgg*. 



Hence the propagation equation can be written as: 



d 



1 = _ i2n ^i N ^. 
z c at* 



(A4) 



In the following discussion it is useful to express the field 
amplitude £j in terms of photon flux r]j, Eq. (|§|), and 
phase ipj. Sep arat ing the real and imaginary parts, we 
find from Eq. ( A4): 



dz 
dz 



dH' 
dtpj ' 

dH' 

drjj ' 



(A5) 



These equations have the form of Hamilton equations of 
classical canonical mechanics with action and angle vari- 
ables r/j, and ipj, "time" z, and the Hamiltonian function 
H' = ^N\ . 

One can see from the eigenvalue equation (|27|) that 
Ao and, hence H', depend on the field phases <pj only 
through the relative phase ip. Therefore, we have: 



dH' _ 2 &H _ dH' 

dipi dip 2 dtp 3 



= 2 



dH' 

dip 



(A6) 



An immediate consequence of this symmetry of H' is the 
existence of constants of motion. Substituting the above 
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equations (A6) into the first line of Eqs. ( A5) yields the 
well-known Manley-Rowe relations 



dm = _ 2 ^m = 2 ^??3 
dz dz dz ' 



(A7) 



which correspond to two independent constants of mo- 
tion: 



Vi + 2??3 = Via + 2? 730, 

m + 2 ??2 = vw + 27720- 



(A8) 



Here 77^0 = Vj( z = 0) are the photon flux values at the 
entrance to the medium. Taking into account the multi- 
photon resonance condition (Bq), one finds furthermore 
that the total intensity of the elm. fields is con served: 
h + 1-2 + h = const(z). The relations Eq. (A8) enable 
us to re-write rjj as: 



r)i(z) = i]io - 2J(z), 
m{z) = V20 + J(z), 
= V3Q + J(z). 



(A9) 



The function J(z) characterizes the amount of energy 
exchange between the waves and has the initial condition 
J(z = 0) = 0. 

Thus the original problem with six amplitude and 
phase variables can be reduced to two variables J and 
ip by a canonical transformation. This leads to 



dJ_ 

dz 
dip 

dz 



&H 

dip ' 
dU 

a/' 



with new Hamiltonian function 



U = -NX + AkJ = -NX. 
2 2 



(A10) 
(AH) 

(A12) 



As can be seen from Eqs. ( A12] ) and (p7f), Ti (or A) does 
not depend on the coordinate z explicitly. Therefore, 
TL (or A) is a fourth constant of motion expressing the 
conservation of the energy density of the medium with 
respect to z. 

To solve the remaining two equations of motion for 
J(z) and ip(z), the Rabi-frequencies Qj are expressed in 
terms of rjjQ and J, and the characteristic equation ( p^ ) 
is written in the form 

G(A, J) = g(J) cos ip. (A13) 

Differentiating both sides with respect to <p yields 

dG _ dGdX 
dip d\ dip 



-— = -gsinip = ±\fg i 



G 2 . 



Substituting this relation into Eq. (A10), we find 



dJ_ 

dz 



± 



2 dG/dX ' 



(A14) 



The choice of the sign in Eq. ( A14) depe nds o n the sign 
of sinip at z = 0. Integration of Eq. ( A14 ) gives an 
implicit solution for J(z): 



± — z 
2 



dG(J') 



dJ' 



dx 



(A15) 



(J') - G* (J') 

Both functions g 2 — G 2 and dG/dX are polynomials in J: 
g = -2/ii^//i 2 ^3 



(A16) 



G 

dG 
~dX 



x \J (vw - uf (V20 + J) (mo + J), 

3 

Go + A ™ jm ' ( A17 ) 



(Ais) 



Therefore, equation ( A14 ) describes a one-dimensional fi- 
nite motion of a pendulum in an external potential. The 
solution is in general given by some combination of ellip- 
tic functions jl5j with parameters determined mainly by 
the roots J n of the polynomial equation: 



g 2 (J) - G 2 (J) = 0. 



(A19) 



The eigenvalue A is a constant of motion (cf. Eq. (A12)), 
and c an thus be found from the characteristic equation 
( A13| ) with parameters taken at the medium entrance 
z = 0: 



G (A) = g(z = 0) cos ip{z = 0). 



(A20) 



Thus, we have reduced the pr opaga tion proble m to s olv- 
ing two algebraic equations: ( A20 ) for A and ( AIS ) for 
the roots J n . If this can be done explicitly, the Hamilto- 
nian method provides an analytical solution to the prop- 
agation problem. But even if an explicit solution is not 
possible, it considerably simplifies numerical calculations. 
The physical meaning of the coefficients A m and a m can 



be drawn by considering the canonical equation (All) for 
the relative phase: 



d<p 
dz 



NdX 
~2dJ 



N dG/dJ _N A 1 + 2A 2 J + 3A 3 J 2 
YdG/dX ~ 



! ao + ai J + a^J 2 

(A21) 

One recognizes that the A m and a m describe the lin- 
ear and nonlinear refraction coefficients of the medium. 
E.g. if J is sufficiently small, t he fi rst term NAi/2ao 
on the right-hand side of Eq. ( A21 ) can be identified 



with the phase mismatch induced by the linear refrac- 
tion, including both contributions from the three-level 
interaction and the residual mismatch Ak. The second 
term (N/2clq) (2A 2 — aiAx/ao) J in the expansion over 
J corresponds to the phase mismatch due to Kerr effect, 
and the next terms are responsible for the higher-order 
contributions. 
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